Exact solution of the anisotropic special transition in the 0(n) model in 2D 



Jerome Dubail 1 , Jesper Lykke Jacobsen 2 ' 1 and Hubert Saleur 1 ' 3 
1 Institut de Physique Theorique, CEA Saclay, 91191 Gif Sur Yvette, France 
2 LPTENS, 24 rue Lhomond, 75231 Pans, France and 
3 Department of Physics, University of Southern California, Los Angeles, CA 90089-0484 

(Dated: September 16, 2009) 

The effect of surface exchange anisotropics is known to play a important role in magnetic critical 
and multicritical behavior at surfaces. We give an exact analysis of this problem in d = 2 for the 
0(n) model by using Coulomb gas, conformal invariance and integrability techniques. We obtain 
the full set of critical exponents at the anisotropic special transition — where the symmetry on the 
boundary is broken down to 0(n\) x 0(n — n\) — as a function of n\. We also obtain the full 
phase diagram and crossover exponents. Crucial in this analysis is a new solution of the boundary 
Yang-Baxter equations for loop models. The appearance of the generalization of Schramm-Loewner 
Evolution SLE K , P is also discussed. 
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The study of boundary critical phenomena is relevant 
to a very large number of physics problems. These in- 
clude, on the condensed matter side, the critical behavior 
of magnets and alloys with free surfaces [l[ , adsorption of 
fluids or polymers on walls and interfaces 0, but also — 
through a series of by now well known mappings — the 
Kondo and other effects in quantum impurity problems 
0, Q]. On the high-energy physics side, apart from the 
old problem of studying field theories on manifolds with 
boundaries, the more recent developments inspired by 
string theory have put boundary effects squarely in the 
limelight. This includes the physics of D-branes 0], and 
of course the celebrated AdS/CFT conjecture [6J. Bound- 
ary effects are also central to many recent developments 
in the study of geometrical critical phenomena and the 
Schramm Loewner Evolution approach [7J. 

Going back to the critical behavior of magnets, surface 
effects can give rise to a bewildering array of physical 
effects [H, especially when combined with finite size ef- 
fects . The experimental activity 1(J II , HJ has been 
steady. 

In interpreting experimental results, a natural question 
concerns the effects of surface anisotropy [la ]. While it 
was quickly understood that such effects are irrelevant 
near the ordi nary transition in the n-vector model 



further study [151 ] showed that they are relevant near spe- 
cial transitions, and that they lead, in high enough di- 
mension, to the emergence of new "anisotropic special" 
transitions. The corresponding family of new boundary 
critical points was extensively studied in d = 4 — e (see 
(l^ l for recent work), and in the 1/n expansion 17]. The 
presence of anisotropy on the boundary is particularly 
interesting in d = 3 since, while only ordinary transi- 
tions are observable for isotropic systems with continuous 
symmetry, the presence of an easy axis on the boundary 
allows for an (anisotropic) special transition, which has 
been studied to second order in e 

We present in this Letter a complete solution of the 



problem in d = 2, in the context of a geometrical re- 
formulation of the n-vector model as a loop model fljj . 
Even though the "surface" is here one-dimensional — and 
so strictly speaking cannot order for integer n — this refor- 
mulation in fact exhibits all the physics of the transitions 
in higher dimension. In particular, we fully recover the 
phase diagram conjectured in Moreover, the loop 

formulation permits to treat n as a real variable, and the 
limits n — > and n — > 1 give access to physical results 
for polymers and the Ising model. 

The case d = 2 is interesting for other reasons as well. 
One concerns the classification of all possible conformal 
boundary conditions for "non-minimal" conformal field 
theories (CFTs) — such as the loop models — with po- 
tential applications to string theory or 2D quantum grav- 
ity [20|. Another has to do with the general program of 
understanding all critical exponents in geometrical mod- 
els, as well as their relations with (variants of) the SLE 
formalism [2l|. 

Loop model. The classical 0(n) model in d dimen- 
sions is defined, initially, by placing on each lattice site 
i a vector spin Si with components S!f for fi = 1, . . . , n. 
Along each link (ij), neighboring spins interact through 
the Boltzmann weight exp(a;S'i • Sj). At high tempera- 
tures (low x) this can be replaced by 1 + xSi ■ Sj . This 
replacement does not modify the long-distance behavior, 
up to and including the critical point x c . The partition 
function then reads 



Z = Trn 

(ij) 



1 



cSi • S- -: 



(i) 



where (ij) denotes the set of nearest neighbors. The 
trace over spin configurations can be normalized so that 
Tr SfSJ = %V 

Expand now the product in |T]), and draw a monomer 
on link (ij) when the term xS^Sj is taken. The trace 
of terms containing an odd power of any S 1 " vanishes 
by the symmetry Si — » — Si. Specializing to a trivalent 



2 



lattice, non-vanishing terms can then only contain 
to the power or 2. The monomers drawn hence form 
configurations C of non-intersecting loops. Summing over 
components then yields 



x x n N 



(2) 



where X (resp. N) is the number of monomers (resp. 
loops) in C. In this formulation, n can be considered a 
continuous variable. When x /* x c , the average length 
of a loop passing through the origin diverges. 

Surface critical behavior. When a boundary is present 
(see Fig. [IJ, boundary sites have fewer neighbors than 
bulk sites, so we can expect boundary spins to be dis- 
ordered at and slightly above x c , where bulk spins start 
ordering. In the loop picture, this means that boundary 
monomers are less probable than bulk monomers, and in 
the continuum limit loops will avoid the boundary. This 
is the ordinary surface transition Ord. 
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FIG. 1: 0(n) loop model with a boundary. Boundary loops 
pass through one or more boundary sites, and can be of type 
1 (drawn blue) or 2 (red). Bulk loops are drawn black. 
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Consider now the model Z = "}2c 
extra weight w each time a loop passes through a bound- 
ary site (see Fig. [T]). When w > 1, loops are attracted 
to the boundary. At a critical value w c , this attraction 
precisely compensates the lower number of neighbors, so 
that bulk and boundary spins order simultaneously. This 
is the special surface transition Sp. All other values of 
w flow towards 1 or oo, the latter being the extraordi- 
nary transition Ex in which a single loop occupies the 
whole boundary. The set of boundary monomers has a 
non-trivial fractal dimension, < cff < d— 1, only at Sp. 

For d = 2 the above references to spin ordering do not 
make sense due to the Mermin- Wagner theorem. The 
transitions Ord and Sp nevertheless exist in the loop 
model. On the honeycomb lattice, for —2 < n < 2 , 

1/2 



x c = (2 + ^/T^n)- 1 ! 2 [22j and w c = 1 + 2/V2^n [2 
Surface anisotropy. Motivated by experiments ll 
and theoretical developments for d > 2 [3, EH, H, 1 



we now allow for anisotropy in the boundary interac- 
tion, breaking the symmetry down to 0{n\) x 0(77,2). 
The effect in the loop model formulation is immedi- 
ate. We call type 1 (drawn blue in Fig. [T]) a loop 
for which /_t = 1, . . . ,m, and type 2 (red) a loop with 
/i = ni + 1, . . . , n. For boundary loops, the sum over \i 
is done separately for each loop type, whereas for bulk 
loops the summation is complete as before. The fugacity 
of boundary loops of type 1 (resp. 2) is thus m (resp. 
?i2 = n — 774). This leads to 



E 



X w 



w 



n n-t 



Tin 



(3) 



where now N is the number of bulk loops only, and we 
have introduced type-dependent weights w\ and wi for 
each time a loop passes through a boundary site. [The 
total weight in Fig. Q] is the product of that of the blue 
(resp. red) boundary loop, x 2A w\ni (resp. x 10 W2ri2), and 
that of the bulk loops (black), x 6 nxx 10 n.] We henceforth 
consider m a continuous variable, with < m < n. 

Integrability. As in previous work 2^, 23|, our exact 
results for the honeycomb loop model are derived from 
the integrability of a related model on the square lattice. 
It is defined by the vertex weights 




where vertices related by horizontal and vertical reflec- 
tion have been drawn only once. This model is integrable 
in the bulk — i.e., the weights satisfy the Yang-Baxter 
equations with spectral parameter u — when [22j 



uji = sin 2i/) sin 3't/' + sin it sin(37/? — u) 

0J2 = sin 2ip sin(3"0 — u) 

^3 = sin 2-0 sin u 

ll>4 = sin7isin(37/; — u) 

W5 = sin(2 , — it) sin(30 — u) 

loq = — sinw sm(ip ~ u) , 



(4) 



and tp parametrizes n — —2 cos 4nj>. Boundary interac- 
tions are introduced as follows 



The following integrable weights — i.e., satisfying the 
boundary Yang-Baxter equation — describe the well- 
known isotropic transitions [23| 



Ord: [3q = sin (|t/> +u) , f3%= (3% = sin (|t/> — u) 
Sp: /3 = cos + u) , $1=02 



cos (|V 



We have generalized this result by finding 
a new anisotropic solution ASi, with n\ = 



3 



sin (4(ji — l)tp) I sin (4ry-0) parametrized by r\: 



fa = sin ((277 
Pi = sin((2?7 
P 2 = sin ((2?7 



u) sin ((2rj 
u) sin ((2?7 
uj sin ((2?7 



(5) 



Another solution AS2 arises from the duality transforma- 
tion that exchanges the two loop types, i.e., rt\ — > n — m. 

When u = -0 we have w 6 = in ((U, and so each vertex 
can be pulled apart horizontally so as to form a pair of 
honeycomb vertices ( < S > ^ < 6 £> ). The weights in ([3]) read 
u>i,2 = \/ Pi,2/{xPa), evaluated 22 



25( in u/2: 



101,2 
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| + ^/l -ni(n-ni) 



(6) 



Continuum limit. The long distance behavior at the 
anisotropic special point AS\ could be inferred by setting 
up the Bethe Ansatz equations corresponding to |5]). An 
alternative and easier route is to use the Coulomb gas 
approach to CFT [Ti|. 

Consider first the bulk theory. Orient each loop inde- 
pendently, and attribute a weight e i7Q /( 27r ' when a loop 
turns an angle a. Summing over orientations, this gives 
n = 2 cos 7. A height field h is defined by viewing ori- 
ented loops as it level lines, across which h — ► h±Tr. This 
is well-known (l8j to renormalize towards a Gaussian free 
field with action 



4ir 



(dh) 



2 d 2 5 



(7) 



The coupling constant g is fixed by requiring that the op- 
erator cos(2/i) conjugate to the discretization of h (and 
hence to the weight n) be strictly marginal. This gives 
g = The central charge c and the Kac critical ex- 

ponents h a fy of primary operators <fr a ,b are parametrized 
by a, b and g: 

o=l- 6 (£^)!, Vt = (»-»',-<«-■>' . (8) 



9 
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It is convenient to define the boundary theory on an 
L x T annulus with period T. Its left rim roles as 
the boundary (cf. Fig. [J). Type 1 loops have weight 
m = sin ((r + 1)7) / sin(r7), with r G (0,tt/7) a new 
parameter. The operators that constrain a boundary 
monomer to be of type 1 or 2 are orthogonal projectors. 
We can write, e.g., the type 1 projector in the basis of 
oriented loops: 



kR 



where the top (resp. bottom) line pertains to the left 
(resp. right) rim of the annulus, and = 2isinr7, 



Jcr = 2isin7. Loops touching only the right rim are 
bulk loops and have a weight n. Requiring weight n\ for 
loops touching both rims fixes ip 1 — ifR = rj. Note that 
our interaction does not conserve the arrow current. 

We assume that the diagonal part induces a flow to- 
wards fixed boundary conditions for the dual field. Wc 
thus end up with a free field with Neumann boundary 
conditions on both rims, and additional weights e ±lri 
for each of the p pairs of oriented half-loops going from 
one rim to the other. This amounts to a height defect 
Ah = 2irp when going around the periodic direction, 
y —* y + T, of the annulus. This can be gauged away 
by writing h(x,y) — h(x,y) + 2irpy/T, where now h is 
periodic. The second term gives rise to 



where we used ([7]) and a Poisson resummation. Integra- 
tion over the first term gives g -1 / 24 / P(q), with modular 
parameter q — e~ 7rT / i and P(q) = Tik>i — 9*0 ■ Using 
now ijHJl gives the exact continuum limit partition func- 
tion in the sector with zero non-contractible loops: 



Z Q (q) 



q -c/24 



nGZ 



,hr., 



(9) 



where the easy limit T ^> L has been used to fix the nor- 
malization. The complete spectrum of critical exponents 
can be read off from ([9]). 

To be precise, ((9]) is valid at the transition AS 2 . This 
can be seen from the limit n± — > n, under which the 
leading exponent h r ^ r in ^ vanishes. This fits with 
AS 2 — > Ord, while ASi — > Sp. However, results for ASi 
can easily be obtained by applying the duality transfor- 
mation to (J9]) . We conclude that the boundary condition 
changing (BCC) operators are <& r , r +i for (ASi/Ord) and 
$ r>r for (AS 2 /Ord). 

Fractal dimensions. A non-contractible loop on the 
annulus is generated by the operator $2,1 0- Confor- 
mal fusing with the BCC operator gives $^,-+1 <8> $2,1 = 
$ r+ i )r +i $,._i /r+ i for AS\. To interpret this, note 
that the first term on the right-hand side is dominant, 
and since w\ > w 2 in ^ this must correspond to the 
insertion of a type 1 (blue) non-contractible loop. The 
second term thus produces a type 2 (red) loop. 

The fractal dimension d^p of the set of type 1 bound- 
ary monomers is conjugate to the operator inserting two 
blue loop strands at the boundary. This is obtained by 
fusing <E> r _|_i jr _|_i with itself, giving the principal contribu- 
tion 3>2r+i 2r+i- Therefore 



/(i) 



1 - h 2r+ i. 2r+ i = 1 - r(r + 1) 



{9 -If 



(10) 



at AS-i. This is non-trivial (0 < dt < 1), just as the 
result [23j df = 1 — ^3,3 = 1 — 2{g — l) 2 /g at the special 



(2) 

transition. We find similarly that dj < at AS\ 
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These results imply the physical interpretation of AS\ : 
type 1 loops are critically attracted towards the bound- 
ary, and type 2 loops retract from it. In other words, 
type 1 (resp. type 2) loops stand at an anisotropic spe- 
cial (resp. ordinary) transition. 

Phase diagram. We are now ready to propose the 
phase diagram of the model ([3]), for x = x c and in the 
regime < n\ < n. See Fig. [2] The fixed points (depicted 
as solid circles) are conformally invariant boundary con- 
ditions, and the double arrows represent flows under the 
boundary renormalization group (RG). 
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FIG. 2: Generic phase diagram for < n\ < n in the rotated 
(wi, W2 ) plane. 

Surface anisotropy must be relevant (resp. irrelevant) 
at the special (resp. ordinary) isotropic transition, since 
the loops see (resp. do not see) the boundary in that case. 
Only in the former case can we expect ri\ to change the 
critical behavior of boundary loops. This agrees physi- 
cally with the conclusions of 14], |l5[ for d > 2. 

More precise evidence for these flows can be found by 
evaluating the boundary entropies Sh = — log g\, for the 
various boundary conditions. By the so-called o-theorem 
24j, 5b increases under the boundary RG flows: the flow 
is from a less stable to a more stable boundary condition. 
gAS! and gAs 2 are related by m — ► n — m. We find [25I ] 



9AS! 



1/4 



sin((r + lh/g) 
sinr7 



sm7 



1/2 



(11) 



,1/2 



to be compared with gord = (2/3) 1 ^ 4 (sin(7/g)/sin7) 
and g Sp = (2/ 5 ) 1 /4 sin (2 7 / ig )/(sin 7 sin(7/ ff )) 1/2 . We 
have thus Sord > Sas!, 2 > ^Sp, as is consistent with 
the flows of Fig. Note that it is possible to flow to 
AS\ by tuning only u>i and not w%, in agreement with 
the interpretation that only type 1 loops stand at AS\. 

Near the point Sp we have w\ — W2 ~ [w c — w) 1 ^ . 
Identifying the operators perturbing in the isotropic and 



anisotropic directions, and using standard scaling argu- 
ments, gives the cross-over exponent </> — (1 — hi ^)/(l — 
/i 3 3). The fact that <f> < 1 for < n < 2 implies the cusp- 
like shape of the phase diagram near Sp. This feature is 
also present 15[ in d > 2. 

Physical realization of SLE Kp . Schramm-Loewner 
Evolutions have a natural generalization with a distin- 
guished boundary point, leading to the two-parameter 
family of processes SLE KtP [7|, |2]J. Its physical relevance 
has however remained unclear so far. 

Our geometric formulation provides a lattice object 
whose scaling limit must be described by such a process. 
Indeed, according to the SLE/CFT correspondance 0], 
type 1 loops (blue) at AS\ correspond to k = i/g and p — 
k (h r +i. r +i — h rt r+i — /i2,i) = ^(4 — k)/2 — 2. The fractal 
dimension of the intersection of the SLE KiP trace with 



the real axis is (flOl): d 



p.) 



(l + p/4) (2 - 8/k - Ap/n). 
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